function J = J_SA( m, param)

    % Calculates J(m) under offer matching
     
    J = (1-param.omega_1).*m./rho_fun(1,'R',param) - param.omega_0./rho_fun(0,'R',param) ...
            + param.J_1*m.^(param.gamma_N_1) + param.J_2*m.^(param.gamma_N_2);

end